Damping by slow relaxing rare earth impurities in Ni8oFe2o 
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Doping Ni8oFe2o by heavy rare earth atoms alters the magnetic relaxation properties of this 
material drastically. We show that this effect can be well explained by the slow relaxing impurity 
mechanism. This process is a consequence of the anisotropy of the on site exchange interaction 
between the 4f magnetic moments and the conduction band. As expected from this model the 
magnitude of the damping effect scales with the anisotropy of the exchange interaction and increases 
by an order of magnitude at low temperatures. In addition our measurements allow us to determine 
the relaxation time of the 4f electrons as a function of temperature. 



The dynamic response of magnetic materials is of fun- 
damental interest as well as essential for various applica- 
tions in modern magnetic storage technology. Often it is 
desirable to tailor the damping parameter and the res- 
onance frequency of magnetic materials independently. 
While the resonance frequency can be controlled rela- 
tively easily by e.g. controlling the saturation magneti- 
zation it is more difficult to change the Gilbert damp- 
ing parameter in a controlled way. Recent experiments 
have demonstrated the ability to modify the damping pa- 
rameter of transition metals and transition metal alloys 
by introducing rare earth (RE) impurities [l|, [2| or 3d- 
and 5d transition metals Q . Unfortunately a convincing 
microscopic understanding of the origin of RE-induced 
damping in metallic alloys is still missing. 

Here, we present experimental results on the magne- 
tization dynamics of thin Ni8oFe2o films doped with the 
lanthanides Gd, Tb, Dy, and Ho. By varying the dopant 
concentration we were able to tune the damping param- 
eter by two orders of magnitude. The dynamic response 
was measured over a wide frequency range (0.5-35 GHz). 
By employing various resonance techniques and configu- 
rations extrinsic and intrinsic relaxation effects are sep- 
arated. This procedure allows us to precisely determine 
the induced damping for the various rare earth dopants 
and, more importantly, to unambiguously identify the 
physical origin if the effect. We find that the slow re- 
laxation of the 4f electron spins of the rare earth atoms 
is responsible for the induced damping. As this mecha- 
nism should also lead to a very strong temperature de- 
pendence of the relaxation we perform temperature de- 
pendent FMR measurements to test the applicability of 
the slow relaxer model. 

A series of 10 and 30 nm thick RE doped Ni8oFe2o 
films was grown by DC magnetron co-sputtering from 
single element targets. A 1 nm thick Ta seed layer was 
first deposited onto the glass substrates. The RE doped 
NigoFe2o films were capped with a 3 nm thick Ta layer 
to prevent oxidation. During deposition an Ar gas pres- 
sure of 2 X 10^"^ mbar was used and the deposition rate 
was about 0.1 nm/s. The film thicknesses of all samples 
were measured by x-ray reflectivity, and the RE concen- 
tration was determined by Rutherford Back Scattering; 
the uncertainty of this method is below 1 at. %. The 



static properties of the samples where investigated by vi- 
brating sample magnetometry and magneto optic Kerr 
effect measurements. In Ni8oFe2o-RE earth intermetal- 
lic alloys the 4f magnetic moments of the RE atoms 
are coupled to the 3d of the Ni8oFe2o moments via the 
intra-atomic 5d orbitals 0, 0- The resulting effective 
5d-3d exchange coupling is antiferromagnetic (AF) lead- 
ing to an antiparallel alignment of RE 4f moments and 
the Ni8oFe2o 3d moments This ferrimagnctic order 

leads to a decreasing saturation magnetization with in- 
creasing RE content. In our samples we observe a linear 
decrease of the saturation magnetization as a function of 
doping for all RE elements with a slope of about 40 Oe 
per atomic percent of RE doping at room temperature 
(RT). All samples have soft magnetic properties with 
small coercive fields (less than 2 Oe) and small uniax- 
ial anisotropy fields (less than 5 Oe). From a structural 
point of view all samples discussed in this letter exhibit 
a poly crystalline fee structure typical for low RE concen- 
trations (below 8%) used here [l|. However clustering 
of the RE atoms does not occur even in the amorphous 
phase at much higher RE concentrations @. 

In ferromagnetic resonance measurements the 
linewidth AH{f) is proportional to the microwave 
frequency / only for Gilbert damping. Two magnon 
scattering due to defects and the superposition of local 
resonance lines due to large scale magnetic inhomo- 
geneities lead to a zero frequency linewidth AH{0) 
If Gilbert damping dominates one has A_ff(0) ^ AH 
and the linewidth at a given frequency can be converted 
into the damping parameter using a = AH^/lo. The 
series of FMR lines shown in Fig. [1] was measured 
at a frequency of 22 GHz using Ho-concentrations of 
— 6 %. This data illustrates the broadening of the 
FMR linewidth as a function of the rare earth concen- 
tration. The pure Ni8oFe2o film exhibits a linewidth of 
approximately 50 Oe whereas AH increases by a factor 
of 20 to 940 Oe for a doping concentration of 6 % Ho. 

The observed linewidth broadening can generally have 
various origins. In order to be able to distinguish contri- 
butions from (i) Gilbert damping, (ii) two magnon scat- 
tering and (iii) sample inhomogeneity we perform FMR 
measurements over a wide frequency range allowing us 
to estimate AH{Q) and thus to verify whether signif- 
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FIG. 1: FMR spectra for four different Ho doping concentra- 
tions measured at room temperature in the in the in-plane 
configuration using a frequency of 22 GHz. The red lines 
show the expected FMR lines given by the saturation mag- 
netization, the g-factor, and the damping constant. In the 
calculation the upshift of the resonance field with increasing 
doping is caused by the decreasing magnetization and increas- 
ing damping constant. The g-factor was determined by out 
of plane FMR measurements and remains nearly unchanged 
up to a doping level of 6 %. Note that for the 6 % sample 
the discrepancy between the expected and the measured line 
position is about 300 Oe. 
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FIG. 2: Frequency dependence of AH for (a) NigoFe2o and Gd 
doped Ni8oFe2o (b) Tb doped Ni8oFe2o (c) Dy doped Ni8oFe2o 
and (d) Ho doped Ni8oFe2o films. All measurements are car- 
ried out in the parallel configuration at RT. 



leant extrinsic contributions are present. Out-of-plane 
FMR measurements further allow one to separate mag- 
netic inhomogeneity (local resonance) and two magnon 
scattering contributions 0,0, as for inhomogeneous sam- 
ples one expects a line broadening in the perpendicu- 
lar configuration compared to the parallel configuration 
(AiJ^ > A77||) For all samples discussed in this 

manuscript the out of plane angular dependence of A_ff 
(not shown) is consistent with Gilbert damping and we 
observe AiJ^ ~ ^^|| • 

Fig. 2 shows the FMR linewidth as a function of fre- 
quency for various concentrations of Gd, Tb, Dy, and 
Ho. The linewidth strongly increases with increasing Tb, 
Dy, and Ho doping concentration, while almost no effect 
is observed for Gd doping (note the 10 x reduced scale 
for Gd doping). For all films at doping levels of 6 % 
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FIG. 3: Damping parameter as a function of the RE con- 
centration at RT for RE=Gd, Tb, Dy, and Ho measured at 
RT. 



and below the linewidth at zero frequency AiJ(O) is very 
small compared to the total linewidth at 22 GHz and 
can be neglected. We conclude that for doping concen- 
trations up to 6 % for all RE dopants the parameter 
a can be easily determined from the slope of the fre- 
quency dependent linewidth. The results are summarized 
in Fig. [21 For Gd doping of 5 % the damping constant of 
the Ni8oFe2o film is not considerably influenced. On the 
other hand Ho, Tb and Dy doped Ni8oFe2o films show a 
very strong dependence of the damping parameter on the 
dopant concentration. With increasing RE concentration 
the damping parameter a increases linearly. From the 
slope of the linear increase we determine the contribu- 
tion to the damping parameter per concentration of RE 
dopant, i.e. a = aNiFo + AaRECRE, where Cre is the RE 
atomic concentration in percent. The values for Aa^E 
are 0.0005, 0.038, 0.036, 0.017 for Gd, Tb, Dy, and Ho. 
We observe that Aq;re for Tb and Dy doping are similar 
and lead to the largest damping, while the value for Ho 
is only about half that of Tb and Dy. For Gd doping 
the contribution is very small and only a consequence of 
the reduced magnetization. This is a striking observation 
considering that earlier experimental 2] and theoretical 
[ll| studies suggested RE induced damping to be pro- 
portional to the orbital moment of the dopants. Based 
on these predictions one should observe the largest effect 
for Ho (L=6) doping, a smaller effect for Dy (L=5) and 
Tb (L=3), and no effect for Gd (L=0). This behavior 
is not observed in our detailed measurements. In Fig. [4] 
the contributed damping parameter Aore is plotted as 
a function of the orbital moment of the dopants. Clearly 
the contributed damping is not proportional to the or- 
bital moment of the dopants. 

In the 1960s Orbach and van Vleck jl2| introduced 
the 'slow relaxing' (SR) impurity model to describe the 
damping in RE-doped Yttrium Iron Garnets (YIG)[l^. 
The essence of the slow relaxer model is the following: 
the 4f multiplet of the RE is split in the moderate ex- 
change field of the 5d electrons. This Zeeman splitting is 
of the order of 10 meV and the levels are hence thermally 
populated at RT. In The anisotropy of the 4f-5d exchange 
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FIG. 4: Comparison of the RE contributed damping Aq:re 
(red bullets) to the anisotropy contribution arising from the 
anisotropic exchange interaction between the 4f moments and 
the conduction electrons if ''''(black line) vs. L [l^. In addi- 
tion values of RE induced damping in Ni8oFe2o are plotted 
from Refs. [§| (gray squares), [ll, H^] (gray star), and [3| 
(gray diamond). In Ref. 01 the damping is given as relaxation 
rates A which were converted into the damping parameter al- 
pha using a = A/(747rM). 



interaction causes a modulation of the 4f exchange spht- 
ting when the 3d moments precess. The thermal popu- 
lation of the 4f levels follows this temporal modulation, 
but is delayed by the RE relaxation time tre. Thermal 
transitions in the 4f multiplet lead to a locally fluctuat- 
ing transverse field acting on the 3d moments via the 
strong 3d-5d coupling if the 4f-5d exchange interaction is 
anisotropic. Using the second fluctuation dissipation the- 
orem one can show that the damping constant is given 
by the Fourier transform of the time correlation of h{£) 
|14| . The time correlation function of the fluctuation field 
can be approximated by {h{t)h{0)) = h"^ exp —t/r lead- 
ing to the following expression for the Gilbert damping 
parameter: 



anE = Cx F{T) 
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l + (u;rRE)2 

where the constant C is given by C — gj^.^^i^^rp 
temperature dependent function F{T) accounts for the 
fact that the precession induced repopulation of the 4f 
levels strongly depends on the temperature Ms is 

the saturation magnetization, Cre is the concentration 
of the RE ions, and A is the anisotropy of the 5d-4f ex- 
change interaction and given by the angular derivatives 
of the 5d-4f exchange energy fifJ 15| . Its magnitude can 
be estimated from the anisotropy contribution arising 
from the anisotropic exchange interaction between the 
4f moments and the conduction electrons as observed in 
metallic rare earth single crystals flol. For a two 
level system one has F{T) = sech^^^ The pop- 

ulation of the 4f levels in RE doped YIG is indeed well 
described as a two level system due to the large crystal 
field spitting. In experiments both ATJre and 5're have 
shown a strong temperature dependence with apeak oc- 
curring at low T when hfl = ksT [H [H [l!]. The 
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FIG. 5: Temperature dependence of the FMR linewidth (AH) 
and the FMR field shift S for a dopant concentration of 2.0 % 
for Tb and Ho, 2.5 % for Dy, and 5 % for Gd. The measure- 
ments where carried out at 10 GHz in the parallel configura- 
tion. The field shift S corresponds to the difference between 
the resonance fields of the doped and undoped samples. Note 
that the data for 2.5 % Dy were multiplied by 4/5 in order 
to be able to compare them to the Tb and Ho 2.0 % data. 
The solid lines represent the expected behavior and were cal- 
culated from Eq. [T] For the temperature dependence tre we 
used Eq. 13 of [l3l. The inset shows the temperature depen- 
dence of Tre as determined from the present FMR measure- 
ments. The red bullets represents earlier measurement of tre 
for Nd doped YIG [13]. 



two level approximation, however, may not be justified 
for RE doped Ni8oFe2o. Due to the absence of a sig- 
nificant crystal field in the metallic alloy 2J-I-1 4f levels 
need to be considered allowing transitions to occur at 2J 
different energies. Therefore considerable broadening of 
the linewidth peak at low temperature (compared to RE 
doped YIG) is expected for RE doped NigoFe2o. 

The real part of EqlT] corresponds to damping and 
causes a linewidth AHre = 5R(q;)^ while the imaginary 
part leads to a negative field shift S're = — 3(q;) — . Pro- 
vided that uj ^ 1/tre the damping is independent of 
the frequency and the resonance field shift is small. The 
negative field shift is a consequence of the time delayed 
damping torque (due to thermal repopulation of the 4f 
levels) leading to an effective longitudinal field. In order 
to predict the relative strength of this effect at a given 
temperature as a function of the RE clement one only 
needs to compare the calculated relative magnitude of 
the anisotropic exchange contribution to the magnetic 
anisotropy [l^ as shown in figlH One finds: = 0, 

K^l = 5.5 X 10^ erg/cm^, = 4.6 x 10^ erg/cm^, and 
= 1.43 X 10* erg/cm'^ [l^. Neglecting the tempera- 
ture dependence in Eq. [1] (it should be roughly the same 
for the various RE elements) one immediately observes 
that the RE induced relaxation should be significantly 
smaller for Ho dopants compared to Tb and Dy. It is also 
apparent from this analysis why doping with Gd with its 
isotropic 4f-5d exchange interaction (S-state) cannot lead 
to additional relaxation. 

The applicability of the SR model for RE doped 
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NigoFe2o can be further tested by verifying whether 
AiJfiE and S're increase with decreasing temperature as 
predicted by Eq[T] Fig. shows the temperature de- 
pendence AiJRE and S're measured at /=10 GHz for 
Ni8oFe2o, Ni8oFe2o:Gd5, Ni8oFe2o:Tb2, Ni8oFe2o:Dy2.5, 
and Ni8oFe2o:Ho2. Indeed for Tb, Dy, and Ho, the data 
can be well described by Eq. [TJ The expected negative 
field shift of the resonance field is clearly observed, see 
Fig. El The temperature dependence of AHre and Sre is 
expected to be similar since it is primarily caused by F(T) 
and TRE(r). Using Eq.[T]the RE relaxation time tre can 
be estimated from the ratio 2|S'RE|/AiJRE = ^tre- From 
the data shown in Fig. [5] for Sre and AiJ^E one can esti- 
mate at room temperature ^ 1 ps and at low tem- 
perature Tre* ^ ~ 3 — 10 ps for Tb, Dy and Ho doping. 
This is in excellent agreement with earlier independent 
measurements [17] for tre observed in RE doped YIG 
as can be seen in the inset of FiglS] The shorter ^ 
causes the field shift to be rather small at RT, cf. FiglT] 
As cijtre ^ 1 for all our measurements Eq. [1] predicts 
a linear frequency dependence of the FMR linewidth in 
agreement with the experimental results shown in Fig. [21 
In addition to the temperature dependence also the 
magnitude of the RE induced linewidth is similar for YIG 
and Ni8oFe2o if one considers the ratio of RE concentra- 
tion to the magnetization. By scaling the YIG results 
from Dillon 19] one expects for a doping level of 2% 
at room temperature a linewidth of a few hundred Oe 
at 10 GHz for RE=Tb and Dy in agreement with the 
present results; see Fig. [SJ It is therefore compelling to 
conclude that the additional damping due to RE doping 
in Ni8oFe2o is caused by the very same slow relaxing im- 
purity mechanism which was originally proposed for RE 
doped YIG. 

Our experimental results sharply contradict earlier 
experimental and theoretical work [llj. Reidy et 
al. performed dynamic measurements using only 



pulsed inductive magnetometry at very low frequencies 
500 MHz). In addition these measurements suffer 
from a large uncertainty with respect to the RE con- 
tent of the samples. The corresponding data points are 
included in Fig. |4l The absolute values for the damp- 
ing constants derived from those measurements are up 
to a factor of 5 lower than our present results. How- 
ever, we would like to point out that our present results 
are in excellent agreement with earlier measurements for 
Tb doped Ni8oFe2o films deposited under similar condi- 
tions as the films used in the present study [l|, [13] (see 
Fig. 2]) . Bailey et al. found a strong dependence of the 
contributed damping on the Ar-pressure during the film 
deposition [l[ , with larger damping observed in films de- 
posited at lower sputter gas pressures. Lower sputter 
gas pressures and lower deposition rates typically lead 
to smoother, more homogeneous films with larger grain 
sizes. Note, that the films used in the present study were 
deposited at even lower sputter gas pressures and depo- 
sition rates than the ones in [l| . The theoretical work by 
Rebel and Hohlfeld 11] is based on orbit-orbit coupling 
effects between RE impurities and itinerant electrons but 
does not consider the slow relaxing impurity model. Re- 
bel et al. solely justify the usage of itinerant electrons 
on the basis of better agreement with the experimental 
data by Reidy et al. [5]. However, their theory results 
in temperature independent Gilbert damping without a 
negative field shift and strictly proportional to the or- 
bital moment of the RE impurities. These predictions are 
clearly at variance with our experimental results. On the 
other hand the present results (temperature, frequency, 
and element dependence of the RE induced damping) can 
be readily explained by the slow relaxing impurity model. 
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